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ABSTRACT 

Context. Shortly after birth, open clusters start dissolving; gradually losing stars into the surrounding star field. The time scale for 
complete disintegration depends both on their initial membership and location within the Galaxy. Open clusters undergoing the 
terminal phase of cluster disruption or open cluster remnants (OCRs) are notoriously difficult to identify. From an observational 
point, a combination of low number statistics and minimal contrast against the general stellar field conspire to turn them into very 
challenging objects. To make the situation even worst, random samples of field stars often display features that may induce to classify 
,— < them erroneously as extremely evolved open clusters. 

Aims. In this paper, we provide a detailed study of the stellar content and kinematics of NGC 6863, a compact group of a few stars 
located in Aquila and described by the Palomar Observatory Sky Survey as a non existent cluster. Nonetheless, this object has been 
recently classified as OCR. The aim of the present work is to either confirm or disprove its OCR status by a detailed star-by-star 
i—{ ' analysis. 

Methods. The analysis is performed using wide-field photometry in the UBVI pass-band, proper motions from the UCAC3 catalogue, 
and high resolution spectroscopy as well as results from extensive iV-body calculations. 

Results. The spectra of the four brightest stars in this field clearly indicate that they are part of different populations. Their radial 
velocities are statistically very different and their spectroscopic parallaxes are inconsistent with them being part of a single, bound 
stellar system. Out of the four stars, only two of them have similar metallicity. The color magnitude diagram for the field of NGC 6863 
does not show any clear signature typical of actual open clusters. Consistently, spatial scan statistics confirms the absence of any 
statistically significant, kinematically supported over-density at the purported location of NGC 6863. 

Conclusions. Our results show that the four brightest stars commonly associated to NGC 6863 form an asterism, a group of non- 
physically associated stars projected together, leading to the conclusion that NGC 6863 is not a real open cluster. 
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from clusters with very different initial populations. Open clus- 
ters originally hosting relatively small numbers of stars (perhaps 
not more than 10 3 ) still arrive to the terminal phase of cluster 
evolution with a significant population of early-type, luminous 
stars. Enhanced preferential retention of massive stars and bina- 
ries is mainly the result of dynamically induced mass segrega- 
tion. In contrast, originally rich clusters (e.g. N > 10 4 ) can sur- 
vive for several Gyr and when they arrive to the twilight of their 
dynamical lives, the brightest stars still in the main sequence are 
late F or early G spectral types. In other words, initially richer 
clusters remain bound for longer time and therefore the main 
sequence tunoff point is fainter by the time that they dissolve. 
This fact makes OCRs from these clusters much more difficult 
to identify than those from smaller clusters. 

Unconnected to Loden's work, a new list of OCR candi- 
dates was produced by Bica et al. (2001); their members be- 
ing selected based on stellar density contrast considerations. 
Unfortunately, subsequent analyses including positional, kine- 
matic, photometric and spectroscopic information failed to con- 
firm even the cluster status of some of the suggested candidates 
(Villanova et al. 2004a; Carraro et al. 2005). Again, negative re- 
sults pointed out an already well known fact, that these objects 
are notoriously difficult to identify. This relatively large frac- 
tion of false positives strongly hints that random samples of field 
stars often display features that may induce to classify them erro- 
neously as extremely evolved open clusters. In spite of their in- 
trinsically challenging nature, the observational search for OCR 
candidates has become a subject of considerable interest as it is 
one of the keys to understand the origin of the field star popula- 
tions of the Galactic disk (see Carraro 2006). 

In this research, we provide a detailed study of the stellar 
content and kinematics of NGC 6863, a compact group of a 
few stars in Aquila described by the Palomar Observatory Sky 
Survey as a non existent cluster. However, this object has been 
recently classified as OCR (Bica et al. 2001; Pavani & Bica 
2007). The aim of the present work is to either confirm or dis- 
prove this result. In this paper, we present new CCD UBVI pho- 
tometry and single epoch spectroscopy, which we combine with 
proper motions from the UCAC3 (Zacharias et al. 2009) cata- 
logue in an attempt to unravel the true nature and dynamical sta- 
tus of this object. This paper is organized as follows: in Sect. 2 
we provide a historical overview of the data previously available 
for this object. Observations, data reduction, and overall results 
are presented in Sect. 3. In Sect. 4 we focus on proper motions. 
The statistical significance of the results is discussed in Sect. 5. 
In Sect. 6 we draw our conclusions. 

2. NGC 6863 in perspective 

This asterism or group of a few stars in Aquila was first identi- 
fied in 1827 by Herschel (1833) as h 2065 using a 18.7 inches 
telescope. The same group of eight Galactic stars was included 
as object GC 4542 by Herschel (1864) in his "Catalogue of 
Nebulae and Clusters of Stars". The "New General Catalogue 
of Nebulae and Clusters of Stars" (Dreyer 1888) first used the 
name NGC 6863, describing the object as a very much com- 
pressed, small cluster including stars from the 19th magnitude 
downwards. No further references to this object are found un- 
til the publication of "The Revised New General Catalogue of 
Nonstellar Astronomical Objects" (Sulentic & Tifft 1973). This 
catalog is a modern, revised, and expanded version of the orig- 
inal NGC. Besides incorporating the many corrections to the 
NGC found over the years, each object was verified on Palomar 
Observatory Sky Survey (POSS) prints and on plates for south- 
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Fig. 1. Palomar Observatory Sky Survey (POSSI) image of the 
area covered by the present study. North is up, East to the left, 
and the image is 12.9 arcmin on a side. This 12 min digitized 
plate was acquired by the Palomar 48-in Schmidt camera using 
the red (O) filter on 1951.50678. The frame displayed has been 
processed using square root scaling. 

ern objects specifically taken for the purpose. In this catalog 
NGC 6863 is described as non existent. The POSSI red (O) 
plate is displayed in Fig. [1] The object is widely regarded as 
an asterism among the amateur astronomical community (see, 
e.g., Allison 2006) but it is listed as a bona fide cluster in 
both the Open Cluster Database Q (WEBDA, Mermilliod & 
Paunzen 2003) and the New Catalogue of Optically Visible Open 
Clusters and Candidates^ (NCOVOCC, Dias et al. 2002). These 
databases are widely used in professional open cluster stud- 
ies. The latest update of WEBDA (November 2009, Paunzen 
& Mermilliod 2009) includes coordinates, a = 20 h 05 m 07 s , 
5 = -03°33'18", / = 38°.278, b = -17°.996, J2000, and di- 
ameter, 2 arcmin, for NGC 6863 but no other data are provided. 
The February 2009 version (v2.10, Dias 2009) of NCOVOCC 
includes NGC 6863 as a well established open cluster with co- 
ordinates, a = 20 h 05 m 07 s , 6 = -03°33'20", I = 38°.2778, 
b = -17°.9962, J2000, and diameter, 3 arcmin, located 1200 pc 
from the Sun and with an age of 3.47 Gyr. SIMBAD0 also classi- 
fies NGC 6863 as cluster of stars and refers to Bica et al. (2001). 
Conversely, the NGC/IC Project0(Erdmann 2009) indicates that 
NGC 6863 is an asterism made of 8 stars. 

An extensive search across published literature on this ob- 
ject reveals only three entries, all of them within the context of 
the Open Cluster Remnant paradigm. Bica et al. (2001) included 
NGC 6863 in their list of dissolving star cluster candidates with 
Galactic coordinates I = 38°.27, b = -17°.99, E(B - V) = 
0.28, and age 3-4 Gyr. Pavani et al. (2007) presented prelimi- 
nary spectroscopy for NGC 6863. Finally and using Two Micron 
All Sky Survey (2MASS, Skrutskie et al. 1997) and UCAC2 

1 http://www.univie.ac.at/webda/ 

2 http://www.astro.iag.usp.br/~wilton/ 

3 http://simbad.u-strasbg.fr/simbad/ 

4 http://www.ngcicproject.org 
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Table 1. Journal of photometric observations in the field of 

NGC 6863 



Field 


Filter 


Time integration 

(sec) 


Seeing 
(") 


NGC 6863 










U 


15, 300 


1.8 




B 


2, 5, 90 


1.9 




V 


1, 3, 30, 30 


1.7 




I 


1,5,30 


1.6 



(Zacharias et al. 2004) data, Pavani & Bica (2007) studied the 
structure, proper motions, and CMD distribution of this object 
to conclude that NGC 6863 is an open cluster remnant of old 
age, 3.5+0.5 Gyr, its centre and characteristic size as in Dias 
(2009), see above. An obvious over-density of bright stars can 
be seen in the central region of Fig. [1] however and following 
(e.g.) Odenkirchen & Soubiran (2002), an overdensity of stars 
does not necessarily imply the existence of a physical ensemble. 

3. Observations and reduction 

Observations were performed in two different runs. The photo- 
metric one was carried out at Asiago (Italy) on December 2001, 
whereas the spectroscopic run was completed at ESO La Silla on 
September 2008. Details of both runs are given in the following 
sections. 

3.1. Photometry 

Observations were carried out with the AFOSC camera at the 
1.82 m Copernico telescope of Cima Ekar (Asiago, Italy), in 
the photometric nights of December 17 and 18, 2001. AFOSC 
samples a 8'. 14 x 8'. 14 field in a IK X IK thinned CCD. The 
typical seeing was between 1.8 and 2.3 arcsec. The details of 
the observations and data reductions are provided in Villanova 
et al. (2004a), where we reported on NGC 5385, Collinder 21 
and NGC 2664, which were observed in the same run. Basic in- 
formation on the NGC 6863 photometric run is given in Table 

m 

From our photometric data we construct color-magnitude di- 
agrams (CMD) for several color combinations. They are shown 
in Fig. [2] Here and with filled symbols we indicate the stars 
defining the asterism at the centre of the NGC 6863 field. These 
four stars produce the appearance of an overdensity in DSS im- 
ages. All of them have UBVI measurements with the exception 
of star #4, which does not have U. The CMDs for the detected 
stars clearly do not show any distinctive feature that can lead us 
to think about a group of stars with common properties, a phys- 
ical ensemble as in the case of, e.g., NGC 1901 discussed in 
Carraro et al. (2007). To guide the eye, and better clarify our con- 
clusion, we over-impose in the same Figure a solar metallicity 
isochrone taken from the Padova database (Girardi et al. 2000), 
and adopting the set of parameters from Pavani et al. (2007): 3.5 
Gyrs for the age, E(B-V) = 0.28, and (m-M) ~ 10.0. 

Our photometric material simply suggests that what we see 
in the direction of NGC 6863 is nothing but random Galactic 
field. 

3.2. Spectroscopy 

The spectra of the four brightest stars in NGC 6863 were col- 
lected at La Silla Observatory on 2008, September 18, with the 



> 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

- 




/ 

/ 

/ 


1 1 


1 1 1 1 1 1 " 1 1 '/" 1 1 " 
/ 

/ 


_ _ 

_ 




/ 

/ 

o ' 
10 




/ 

— / — 

o / 

/ ° 


<" " o 








- \ V ° 






/ 

1 > • 

*• o 




- \ 


- ! o ? 1° 




\ S 8 




\ o° S 






<s> 




, Oo 






•V 

6 




— ' O 
i° 8 ° 










\ o 


p k? °° ° 

\ 








Jfr ^ 


o \ 








: ° Jte : 


\ 

_ \ _ 








- ^ \ - 


\ 

\ 








- o V o°\ - 


\ 

\ 




° 1 




\ " 


1 1 1 1 I r i r T I T 1 T T I T 1 1 1 1 




1 




V 



0.5 1 1.5 20 0.5 1 1.5 0.5 1 1.5 2 

(U-B) (B-V) (V-I) 



Fig. 2. CMDs in various color combinations for all the stars hav- 
ing UBVI photometry in the field of NGC 6863. Filled simbols 
indicate the 4 stars which identify the over-density (see Fig. [TJ. 
Super-imposed are a solar metallicity isochrone from Girardi et 
al. (2000) adopting the set of parameters derived by Pavani et al. 
(2007) for this asterism. 

HARPS fibre-fed spectrograph at the Cassegrain focus of the 
3.6m telescope. Exposure times varied between 1800s and 3600s 
depending on target brightness, and the resulting signal-to-noise 
ratio for each spectrum is given in column 2 of Table [2] The 
light contribution from sky background was not negligible for 
our faint targets, and we allocated the second fiber to the sky 
instead of the simultaneous wavelength-calibration lamp, thus 
relying for calibration only on lamp frames acquired during day- 
time operations. However, the very high spectral resolution and 
instrumental stability guaranteed an accuracy still much higher 
than what was required for our aims. 

We first attempted a manual data reduction by means of stan- 
dard IRAf0routines, comprising all steps from bias and flat-field 
corrections to spectra extraction and wavelength calibration. The 
inspection of the results revealed that there was no improvement, 
in terms of S/N and spectral quality, with respect to the output of 
the online pipeline at the telescope. The wavelength calibration 
was also equivalent within calibration errors, and no systematic 
effect was detected. Hence, we adopted the pipeline-reduced, 
extracted and calibrated spectra. The extracted sky spectrum 
from the second fiber was subtracted to each target. The order 
#46 of the second fiber falls between the two CCDs, and the 
corresponding science order was also suppressed, generating a 
gap between 5258 and 5336 A. Spectra were convolved with 
a Gaussian filter to decrease the resolution to about R=38 000, 
then re-binned in blocks of 4 pixels and re-sampled to a constant 
wavelength step of 0.01 A. We thus increased the S/N by a fac- 
tor of two, while still keeping the resolution required for precise 

5 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. 3. A portion of the collected spectra, centred on the strong 
Nal doublet. An integer constant have been added to shift spectra 
vertically. 



radial velocity and abundances measurements. Resulting spectra 
were divided order-by-order by the blaze function obtained by 
the pipeline from fiat field frames, and any residual curvature of 
orders was eliminated fitting the continuum with a second-order 
polynomial. The orders thus normalized were finally merged to 
produce the final spectra. A section of the final spectra across 
the strong Nal doublet (5890-5893 A) is presented in Fig. [3] The 
whole useful spectral range goes from 4000 to 6800A. 

3.2.1. Radial velocities 

Radial velocity (RV) measurements were performed with a 
cross-correlation (CC) technique (Tonry & Davis 1979) using 
the fxcor IRAF task. As template we used a synthetic spec- 
trum of solar metallicity drawn from the library of Coelho et al. 
(2005), with temperature and surface gravity intermediate be- 
tween the values expected for program stars (4750 K and 2.0 
dex, respectively). We verified through repeated measurements 
that different assumptions on template parameters in the range 
T e ff=4000-5500 K and log §=1.0-3.0 caused negligible changes 
in the resulting velocities (<0.2 km s ). The whole useful spec- 
tral range was used in the CC, excluding the gap at 5300 K and 
the head of the telluric band at the red end of the range. In order 
to fix the zero-point of the measurements, this band was cross- 
correlated with a twilight solar spectrum acquired during a dif- 
ferent run. Within 0.1 km s the CC returned the solar RV at 
the dates of observations^, indicating that no correction for sys- 
tematic effects was required. Finally, observed RVs were cor- 
rected to heliocentric velocities. Errors have been estimated as 
the quadratic sum of the most relevant sources, i.e. CC error 
(0.2-0.4 km s ), the uncertainty in zero-point definition (0. 1 km 



s ), and choice of the synthetic template (0.2 km s '). Results 
are given in column 7 of Table [2] 

3.2.2. Atmospheric parameters 

Effective temperature (T e ff), surface gravity (logg), and metal- 
licity of stars #4 and #5 were measured with the ionization and 
excitation balance method. In brief, equivalent widths (EWs) of 
a great number of Fel and Fell lines are measured, then T e ft 
and microturbulent velocity (v,) are determined imposing that 
Fel abundances (e(Fel)) of each individual line show no trend 
with excitation potential (EP) and EW, respectively, while log g 
is fixed by the requirement that Fel and Fell lines return the same 
iron abundance. We adopted the general iron line list of Moni 
Bidin (2009), which comprises atomic parameters for about 180 
Fel and 40 Fell lines, with solar oscillator strengths and a line- 
by-line definition of the correction term to the Unsold (1955) 
approximation of the collisional line broadening. However, we 
rejected the measurement when the line fit was not satisfactory 
or the stellar continuum was not easily defined, and only ~70 
iron lines were finally used for each star. EWs were measured 
with a Gaussian fit with a dedicated IRAF script. Iron abun- 
dances were calculated with the LTE code MOOCfl (Sneden 
1973), with model atmospheres obtained interpolating from the 
Kurucz (1992) grid. The results are given in Table [2] We adopt 
the error estimates of Moni Bidin (2009), who compared the re- 
sults with literature values for thirteen stars with multiple previ- 
ous measurements. However, we improved the problematic fit of 
Fell lines, found at shorter wavelengths and affected by heavy 
line crowding, by simultaneously fitting all features in a region 
centred on the line of interest. Measurements on some test stars 
indicate that the scatter of log g is thus reduced, and we adopt 
here a corresponding lower error. 

EW measurement for stars #7 and #8 was possible for too 
few lines (<30) and with large uncertainties, because of the 
low S/N. Consequently, the atmospheric parameter determina- 
tion failed. For example, for star #7 there was no physically 
meaningful combination of parameters that permitted to zero the 
trend of e(Fel) with EWs and EP simultaneously. We therefore 
elaborated a different strategy. The wings of Balmer lines are 
a good indicator of temperature in the range T e ff =5000-6500 K 
(Fuhrmann et al. 1994), and an explorative inspection of syn- 
thetic spectra revealed that they are quite insensitive to log g. We 
fitted the wings of H a on target spectra with synthetic templates 
at different temperatures in steps of 250 K, and results were fit- 
ted using a second-order polynomial to determine the minimum 
of the ^ 2 -T e ff relation. The estimate was repeated within a grid 
of values of metallicity and gravity (logg from 1 to 4, [Fe/H]=0 
and -0.5) to check the extent of systematics, but only small ran- 
dom differences were found. We will assume the mean of these 
repeated measurements as the final result, and their rms as the 
associated error. We first tested the procedure on star #1, obtain- 
ing 5400±65 K. The excellent agreement with our estimate ob- 
tained with the EW method indicates that the results are reliable. 
It was not possible to repeat the test on star #4, because Balmer 
lines rapidly vanish at decreasing temperature, and their reduced 
wings become insensitive to temperature changes. In fact, for 
this cool star the^- 2 -T e ff was flat below about 4300 K. 

A similar routine was prepared to estimate logg, fitting 
the wings of three gravity-sensitive lines: the Cal line at 6162 
A (Edvardsson 1988, Katz et al. 2003), the redder line of the 
Mglb triplet at 5182 A(Kuijken et al. 1989), and the Nal doublet 



6 http://eclipse.gsfc.nasa.gov/TYPE/TYPE.html 



7 freely distributed by C. Sneden, University of Texas, Austin 
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Table 2. Stellar data derived with spectroscopic analysis 



ID 


s 

N 






l°gg 


r— i 

L H J 


v, 


RV 


E(B - V) 


d 






K 




dex 


dex 


km s~' 


km s" 1 


mag 


pc 


4 


25 


4110+ 


65 


1.05±0.20 


-0.62+0.08 


1.28+0.20 


-55.1+0.3 


0.23+0.05 


6800+1100 


5 


15 


5440+ 


65 


3.55±0.20 


-0.09±0.08 


0.65±0.20 


21.6±0.5 


0.30±0.05 


1270±280 


7 


10 


6100+ 


75 


4.45±0.20 


-0.06±0.09 


1.32+0.15 


35.3±0.4 


0.34±0.04 


560±50 


8 


10 


5320+ 


65 


3.85±0.10 


-0.41±0.10 


1.24+0.15 


-10.6+0.5 


0.26+0.05 


740±90 




5180 5181 5182 5183 5184 5185 5186 5187 



I 

Fig. 4. The Mgl line of stars #7 and #8, with overplotted (solid 
lines) synthetic profiles of templates with corresponding tem- 
peratures. The plotted models are for log g=3.0 (narrower, upper 
profiles), log g=4.0 (middle lines), and log g=5.0 (wider, lower 
profiles). The spectra of star #8 was vertically shifted adding a 
unit constant. 



at 5890-5893 A . The temperature was fixed to the previously- 
determined value, interpolating the spectra using the synthetic 
library of Coelho et al. (2005), while both solar metallicity and 
[Fe/H]=-0.5 were used. The three lines were fitted indepen- 
dently, resulting in six gravity measurements per star. For stars 
#4 and #5 we obtained log g= 1.00+0. 10 and 3.70+0.10 respec- 
tively, confirming again the reliability of the results. In Figure|4] 
the Mgl line of stars #7 and #8 is shown. As it can be seen, syn- 
thetic line profiles change noticeably with gravity, and a rough 
estimate of stellar gravity can even be done by visual inspection. 

In order to fix the microturbulence velocity, we considered 
the equations available in the literature that estimate this param- 
eter as a function of T e ff and log g. For star #7 our choice fell 
on the relation presented by Reddy et al. (2003), that was cal- 
ibrated using stars of the same spectral type as the target and 
presents a small quoted error (0.15 km s~'). Star #8 is cooler 
than the temperature range where this equation was calibrated. 
The relation proposed by Gratton et al. (1996) returns a much 
lower v, for cool stars, but it was calibrated on a mix of dwarfs 
and giants in a large range of temperatures, and the quoted un- 
certainty is much larger (0.30 km s~'). As no approach was fully 
satisfactory because both equations had a characteristic weak- 



ness, we decided to adopt the mean value as the best estimate of 
v, for star #8. However, different assumptions have only limited 
effects, because both relations return a metallicity within lcr of 
the final estimate. 

The metallicity of stars #7 and #8 was measured from the 
EWs of the few iron lines for which a reliable fit was possible. 
The number of available lines was small (about 15), but their fit 
was good and the atmospheric parameters had already been de- 
termined. The rms of [Fe/H] from single lines was about 0.15 
dex, resulting in a formal error on the mean of about 0.04 dex. 
This was quadratically summed to the variations of [Fe/H] in- 
duced by an increase/decrease of one parameter (T e ff, log g, v,) 
by lcr. 

3.2.3. Reddening and distances 

We finally compared the position of program stars in the T e ff- 
log g plane with Yale-Yonsei isochrones (Yi et al. 2003) of cor- 
rect metallicity. The star position in this plane does not depend 
on reddening or distance, but it can depend on age. For each 
star we adopted the isochrone that matched the stellar param- 
eters. However, this gives a reliable age estimate only for sub- 
giant stars (#5 and #8), while at the tip of the RGB (star #4) 
and in the MS (#7) age is ill-defined because isochrones are too 
closely spaced. For stars #5 and #8 we obtained an age estimate 
of 2.3± 1 .0 and 12+3 Gyr respectively, where the errors were de- 
fined by the most extreme isochrones compatible with the la- 
errors in temperature and gravity. We calculated the color ex- 
cess E(V - I) by comparison of observed and theoretical colors. 
E(V - I) was then translated into E(B - V) using the relations 
found by Cardelli et al. (1989). The two estimates of redden- 
ing usually differed because of observational errors, hence we 
adopted as final value their weighted mean with associated error. 
Results are presented in column 8 of Table [2] Uncertainties on 
color excesses were calculated from propagation of errors on ob- 
served and intrinsic colors, the latter being estimated as the vari- 
ation induced by a change of T e ff by + 1 cr. In fact, we verified that 
the intrinsic color is minimally affected by the errors on the other 
parameters (logg, [Fe/H] and age). All but one of the values are 
very similar, and the mean reddening is E(B - V)=0.28. The ex- 
tinction map of Schlegel et al. (1998) returns E(B - V)=0.295 
in the direction of NGC 6863, in excellent agreement with our 
estimate. This also confirms that the temperatures adopted for 
program stars are substantially correct. 

From isochrones we also estimated the absolute magni- 
tude of each object, and we calculated the distance assuming 
E(B - V)=0.28 and A v =3. 1 E(B - V). The error on M v was as- 
sumed as the quadratic sum of the differences induced by the 
variation of each parameter (T e ff, logg, [Fe/H] and age) by lcr, 
and the uncertainty on distance was calculated from propaga- 
tion of errors. Results are presented in column 9 of Table|2] It is 
worth noticing that the distance of the MS star #7 is well-defined, 
while errors on surface gravity induce large uncertainties on the 
distance of subgiants (#5 and #8). 
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Fig. 5. Extended proper motion analysis. Top panel: stars having 
12 < K < 14. Middle panel: stars having 10 < K < 12. Bottom 
panel: stars brighter than K = 10. 
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Fig. 6. Proper motion analysis. UCAC3 stars within 3 arcmin 
from the centre of the object. 



there is no isochrone with [Fe/H]=-0.08 and E(B - V)=0.28 
that fits the stellar position of both stars in the CMD. This differ- 
ence excludes the possibility that they are part of the same stellar 
system. As a conclusion, the spectroscopic analysis reveals that 
the four brightest stars of NGC 6863 are not physically related; 
they are just a random alignment along the line of sight. In other 
words, they are an asterism. 



It is interesting to note that our distance and age for star #5 
agree well with the values cited by Dias (2009) for NGC 6863 
(1200 pc, 3.47 Gyr), but this is true for this star only. 

3.2.4. Spectroscopic results: discussion 

Metallicity, distance and radial velocity of star #4 clearly show 
that it is a distant background star without any connection with 
the others. At a height z=2.1 kpc from the Galactic plane and 
with intermediate metallicity, it appears as a typical thick disk 
red giant. To a lesser extent, a similar conclusion can be drawn 
for star #8 also. Its properties are less extreme, but its RV is 
very different from those of the two remaining stars, and its 
metallicity is lower by more than 3cr. Stars #5 and #7, on the 
other hand, have very similar metallicity. Their RVs differ by 
about 14 km s _1 , but this fact alone does not exclude that the 
two objects can be part of a stellar system, because they were 
observed in a single epoch only, and one or both could be bi- 
naries. It is worth pointing out that indeed a very high binary 
fraction is expected in OCRs (de la Fuente Marcos 1997, 1998), 
and among the 13 confirmed members of NGC 1901, to date 
the only spectroscopically-confirmed OCR, Carraro et al. (2007) 
measured at any epoch an internal RV dispersion between 16 and 
20 km s _1 . However, the spectroscopic distance of stars #5 and 
#7 still differ by 2.5<x. This is also consistent with the fact that 



4. Proper motions 

Following Carraro et al. (2007), we extracted the proper mo- 
tion components in a 20-squared arcmin field around NGC 6863 
from the UCAC3 catalogue (Zacharias et al. 2009), and con- 
structed 3 vector-point diagrams as a function of the K magni- 
tude (UCAC3 is cross-correlated with 2MASS). The results are 
shown in Fig. [5] where the bottom panel is restricted to stars hav- 
ing K < 10, the middle panel to stars in the range 10 < K < 12, 
and the top panel to stars having 12 < K < 14. Stars in the 
lower panel appear to clump around p a = -2.0 [mas/yr], pg = - 
7.5 [mas/yr]. Identifiable, consistent clumps are also observed in 
the other two panels. Is this apparent clump the kinematic signa- 
ture of NGC 6863? 

Following Dias (2009) and restricting our analysis to the cen- 
tral region, 3 arcmin around the centre of the cluster, the results 
are rather different, see Fig. [6] Clumping, if any, is observed 
around p a = -1.5 [mas/yr], p$ = 1.5 [mas/yr]. This diagram in- 
cludes all the stars brighter than K = 14 with UCAC3 data within 
3 arcmin of the centre of the object. In principle, we are using 
data similar (they use UCAC2 not UCAC3) to those in Pavani 
& Bica (2007) but a direct comparison is not possible as they 
do not provide a diagram to illustrate the modulus of the proper 
motion distribution for NGC 6863. They only indicate that its 
characteristic median velocity is the highest in their open cluster 
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remnant sample, 137 km/s. In any case, preliminary visual in- 
spection of Fig. [6] appears to confirm that the stellar population 
located around the commonly accepted centre of NGC 6863 is 
kinematically heterogeneous which is incompatible with open 
cluster status. In the following section, this conclusion is ex- 
plored by using a more detailed statistical analysis. 

As for the four brightest stars discussed above, individ- 
ual proper motions are shown in Table [3] Calculating the 
Galactic space velocity for these stars as described in Johnson 
& Soderblom (1987) with UCAC3 data in the Heliocentric ref- 
erence frame velocity UVW, we obtain results fully consistent 
with our previous spectroscopic analysis. For star #4, U = 92 
km s _1 , V = -320km s , and W = -211 km s _1 , which is con- 
sistent with a halo star in a retrograde orbit. Star #5 has U = 52 
km s _1 , V — -14 km s _1 , and W = 28 km s _1 . Similar numbers 
are found for star #7 with U - 51 km s _1 , V = 2 km s _1 , and 
W = 11 km s _1 . Star #8 is unusual and it could be a thick disk 
member with U — 25 km s _1 , V = 58 km s _1 , and W = 206 km 
s- 1 . 



5. An asterism in Aquila: statistical analysis 

At this point, one may want to know how significant, statistically 
speaking, is the kinematic evidence against NGC 6863 being a 
real open cluster. It is true that the four brightest stars in the 
field of NGC 6863 do not appear to form a physical system but 
what if just one of them is the brightest member of a real but 
poorly populated open cluster like NGC 1901? We are facing a 
spatial data mining problem where the goal is detection of a hy- 
pothetical over-density of kinematic nature (i.e. a small group of 
stars with similar proper motions) but in addition, statistical test- 
ing must be performed in order to determine whether the over- 
dense regions are significant. In order to distinguish the clumps 
that are significant from those that are likely to have occurred by 
chance we use Kulldorff 's spatial scan statistic (Dk defined be- 
low, Kulldorff 1997). In the following, we assume data has been 
aggregated to an N x N grid of squares, where each grid cell is 
associated with a count cy (number of stars with certain proper- 
ties, proper motions in this case) and an underlying population 
Pij, with i, J=l, N. Our goal is to find over-densities: spatial 
regions where the counts are significantly higher than expected, 
given the underlying population. Our objective is to search over 
all square regions and find the region(s) with the highest density 
according to a density measure. Kulldorff 's statistic assumes that 
counts Cy are generated by an inhomogeneous Poisson process 
with mean q py, where q is the underlying expected value of the 
count to population ratio. We then calculate the logarithm of the 
likelihood ratio of two possibilities: that q is higher in the region 
than outside the region or that it is identical inside and outside 
the region. For a cell with count cy and population py, in a grid 
with total count C and population P, we can calculate: 

D Kij = cy log + (C - Ci j) log ^f^- - C log § ( 1 ) 

if dj I pij > C/P, and D Kij = otherwise. Kulldorff (1997) 
proved that the spatial scan statistic is individually more pow- 
erful for finding a single significant region of unusually high 
density than any other test statistic. Unfortunately, performance 
of scan statistics is affected by tessellation geometry (variation 
of cell size and cell shape) as well as aggregation level (average 
cell size). A single actual cluster that is bisected by a large cell 
may appear either as not significant or as two distinct clusters 
depending on whether the large cell is included in the candidate 



zone (Kulldorff, Tango, & Park 2003). In the following analysis, 
special care has been taken to minimize tessellation effects by 
using multiple grid sizes and checking for consistency. On the 
other hand, we have calculated D^j for each cell in the grid, the 
standard deviation, and plot {DKij- < Dk >)/o~ using colour 
maps. To test the reliability of our statistical approach we first 
apply it to a well established OCR, NGC 1901, and then to our 
current object, NGC 6863. 

In order to carry out this analysis we use data from 2MASS. 
The 2MASS All Sky Catalog of Point Sources (Cutri et al. 2003) 
includes data for both NGC 1901 and NGC 6863 and 20x20 
arcmin 2 areas around the coordinates of the accepted centre of 
these objects were extracted. Following Pavani & Bica (2007) 
we use E(B - V) = 3.03xE(J - H), Aj/A v = 0.282 and Aj = 
2.65x£(7 - H) (Rieke & Lebofsky 1985). 



5.1. NGC 1901 

Using UCAC2 data, Carrara et al. (2007) concluded that the 
average proper motions of this genuine OCR were (p a ,p s ) = 
(1.7+1.3, 12.3+2.9) mas/yr. Repeating the analysis with UCAC3 
data for the suspected NGC 1901 members compiled in Table 2 
of that paper we obtain (p a ,p 6 ) — (1.3 ± 1.5, 1 1.9 + 1.5) mas/yr. 
This translates into U = -22.4 ± 8.4 km s _1 , V = -3.0+1.7 
km s _1 , and W = -1.9 + 2.8 km s , which is consistent with 
other recent results (see Bovy et al. 2009, Fig. 21). If we ap- 
ply the spatial scan statistics to the corresponding UCAC3 sam- 
ple (20x20-arcmin 2 around the cluster centre) looking for over- 
densities within 2cr of the average proper motions, we obtain 
the plot in Fig. (right panel). The effective size of the OCR is 
~3 arcmin (Carrara et al. 2007); this is the grid cell size used 
in the analysis. A statistically significant over-density, ~3<x, ap- 
pears at the expected position. The 2MASS CMD is also pro- 
vided, E(B — V) = 0.04 was used in the reduction process. The 
CMD in Fig.|7]is analogous to the one in Pavani & Bica (2007) 
even if a slightly different E(B - V) was used; see their Fig. 7, 
central panel. Photometry does show a defined main sequence 
and the CMD hints of all the characteristic properties of OCRs 
as predicted by /Y-body simulations (de la Fuente Marcos 1998): 
scarcity of low-mass members, high binary fraction (a double 
main sequence is not unlikely), and overall sparse population. 
These results together with the kinematically homogeneous sig- 
nature confirm that the statistical analysis described above gen- 
erates robust results. NGC 1901 is a true representative of the 
terminal stage in the life of an open cluster, the OCR phase. Its 
current properties are compatible with NGC 1901 being what re- 
mains of a relatively small cluster with an original population of 
500-750 stars (Villanova et al. 2004b, Carrara et al. 2007). 



5.2. The NGC 6863 asterism 

If we repeat the same approach looking for over-densities associ- 
ated to the four brightest stars in the field of NGC 6863 we obtain 
Fig- IS Here, the effective size of the hypothetical OCR is ~1.5 
arcmin (Pavani & Bica 2007) and this is the grid cell size used 
in the analysis. E(B -V) = 0.28 was used to plot the CMDs. No 
obvious over-density can be found at the location of NGC 6863. 
It is, therefore, reasonable to conclude that the actual presence 
of an open cluster, remnant or not, at the location of NGC 6863 
is very unlikely. The application of the Kulldorf statistics shows 
that NGC 1901 is a real OCR, whereas it confirms that there is 
no evidence to suggest that NGC 6863 really exists. 
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Table 3. Proper motion data for the brightest stars 



ID 


a(72000) 


(5(72000) 


fi a UCAC3 


Us UCAC3 






(° ' ") 


[mas/yr] 


[mas/yr] 


4 


20:05:08.50 


-03:33:22.7 


2.2±3.9 


-11.9+6.8 


5 


20:05:07.67 


-03:32:59.3 


-8.6±3.4 


-3.6±4.1 


7 


20:05:06.07 


-03:33:29.3 


-12.9±15.2 


-6.1±5.4 


8 


20:05:06.78 


-03:33:26.6 


-49.9±6.1 


35.5±6.0 



6 




Fig. 7. NGC 1901 . Left panel: 2MASS CMD for those stars with proper motions in UCAC3 within 2cr of (ji a , p s ) = (1.3 ± 1.5, 1 1.9 ± 
1.5) mas/yr. Right panel: Kulldorff's spatial scan statistics for the region around NGC 1901 (20x20 arcmin 2 field, UCAC3 data). 
There is an obvious, statistically significant (~3cr) over-density of stars with similar proper motions centred at the accepted coordi- 
nates for this OCR. 



6. Conclusions 

In this paper we have presented detailed photometric, spectro- 
scopic and kinematic evidence against the classification of NGC 
6863 as open cluster. Our results show that the four brightest 
stars commonly associated to NGC 6863 form an asterism, a 
group of non-physically associated stars projected together. The 
spectra of the four brightest stars in this field clearly indicate 
that they are part of different populations. Their radial velocities 
are statistically very different and their spectroscopic parallaxes 
are inconsistent with them being part of a single, bound stel- 
lar system. Out of the four stars, only two of them have similar 
metallicity. As for the underlying open cluster remnant described 
in Pavani & Bica (2007), our analysis of the existing data has 
demonstrated convincingly that there is no base whatsoever to 
consider NGC 6863 as a possible OCR. The color magnitude di- 



agrams for the field of NGC 6863 do not even show any clear 
signature typical of actual open clusters. Spatial scan statistics 
using UCAC3 data strongly suggests that no statistically signifi- 
cant, kinematically supported over-density is present at the pur- 
ported location of NGC 6863 although the proper motion errors 
are rather large. This together with the other issues discussed 
above force us to conclude that the presence of a star cluster of 
any kind associated to the coordinates of NGC 6863 is highly 
unlikely. 

Acknowledgements. CMB acknowledges the Chilean Centro de Excelencia en 
Astrofi'sica y Tecnologfas Afines (CATA), and Francesco Mauro for very use- 
full discussions. GC deeply acknowledges the entire Asiago technical staff for 
the kind night assistance over the whole duration of this project. In prepara- 
tion of this paper, we made use of the NASA Astrophysics Data System and 
the ASTRO-PH e-print server. This research has made use of the WEBDA 
database operated at the Institute of Astronomy of the University of Vienna, 



Moni Bidin et al.: The NGC 6863 asterism 



9 




10 



Moni Bidin et al.: The NGC 6863 asterism 



Austria. This work also made extensive use of the SIMBAD database, op- 
erated at the CDS, Strasbourg, France. This publication makes use of data 
products from the Two Micron All Sky Survey, which is a joint project of 
the University of Massachusetts and the Infrared Processing and Analysis 
Center/California Institute of Technology, funded by the National Aeronautics 
and Space Administration and the National Science Foundation. 



References 

Aarseth, SJ. 1968, Bull. Astron., Ser. 3, 3, 105 

Allison, M. 2006, Star Clusters and How to Observe Them, Springer- Verlag 

London Limited, p. 45 
van Albada, T.S. 1968, Bull. Astron. Inst. Neth., 19, 479 

Bica, E., Santiago, B. X., Dutra, C. M., Dottori, H., de Oliveira, M. R., & Pavani, 

D. 2001, A&A, 366, 827 
Bovy, J., Hogg, D. W., Roweis, S. T. 2009, ApJ, 700, 1794 
Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 245 
Carraro, G. 2006, BASI, 34, 153 

Carraro, G, Dinescu, D. I., Girard, T. M., & van Altena, W. F. 2005, A&A, 433, 
143 

Carraro, G, de La Fuente Marcos, R., Villanova, S., Moni Bidin, C, de La 
Fuente Marcos, C, Baumgardt, H., & Solivella, G. 2007, A&A, 466, 931 

Coelho, P., Barbuy, B., Melendez, J., Schiavon, R. P., & Castilho, B. V. 2005, 
A&A, 443, 735 

Cutri, R., Skrutskie, M. F, Van Dyk, S., et al. 2003, The 2MASS All-Sky Catalog 
of Point Sources, University of Massachusetts and IPAC/California Institute 
of Technology. 

Dias, W. S. 2009, Catalog of Optically Visible Open Clusters and Candidates, v. 

2.10., http : //www . astro . iag . usp . br/~wilton/ 
Dias, W. S., Alessi, B. S., Moitinho, A., & Lepine, J. R. D. 2002, A&A, 389, 871 
Dreyer. J. L. E. 1888, Mem. R. Astron. Soc, 49, 1 (reprinted 1962, London: 

Royal Astronomical Society) 
Erdmann, R. E. 2009, The NGC/IC Project, [http : //ww w . ngcicp roject . org| 
de la Fuente Marcos, R. 1997, A&A, 322, 764 
de la Fuente Marcos, R. 1998, A&A, 333, L27 
Edvardsson, B. 1988, A&A, 190, 148 

Fuhrmann, K., Axer, M., & Gehren, T. 1994, A&A, 285, 585 
Girardi, L., Bertelh, G, Bressan, A., Chiosi, C, 2000, A&AS, 141, 371 
Gratton, R. G, Carretta, E., & Castelli F. 1996, A&A, 314, 191 
Herschel, J. F. W. 1833, Phil. Tran. R. Soc, 123, 359 
Herschel, J. F. W. 1864, Phil. Tran. R. Soc, 154, 1 
von Hoerner, S. 1960, Z. Astrophys., 50, 184 
von Hoerner, S. 1963, Z. Astrophys., 57, 47 
Johnson, D. R. H., & Soderblom, D. R. 1987, AJ, 93, 864 
Katz, D., Farata, F, Aigrain, S., & Micela, G. 2003, A&A, 397, 747 
Kuijken, K., & Gilmore, G. 1989, MNRAS, 239, 605 
Kulldorff, M. 1997, Commun. Statist.: Theory Meth., 26, 1481 
Kulldorff, M., Tango, T., & Park, P. 2003, Computational Statistics and Data 
Analysis, 42, 665 

Kurucz, R. L. 1992, in IAU Symp. 149, The Stellar Population of Galaxies, ed. 

B. Barbuy & A. Renzini (Dordrecht Reidel), 225 
Lada, C. J., & Lada, E. A. 2003, ARA&A, 41, 57 
Loden, L. O. 1987, Ir. Astron. J., 18, 95 
Loden, L. O. 1988, Ap&SS, 142, 177 
Loden, L. O. 1993, Ap&SS, 199, 165 
Mermilliod, J.-C, & Paunzen, E. 2003, A&A, 410, 511 
Moni Bidin, C. 2009, PhD Thesis, Universidad de Chile, 1 
Odenkirchen, M., & Soubiran, C. 2002, A&A, 383, 163 

Paunzen, E., & Mermilliod, J.-C. 2009, Open Cluster Database, 

http: //www.univie.ac . at/webda 
Pavani, D. B., & Bica, E. 2007, A&A, 468, 139 

Pavani, D. B., Bica, E., Ahumada, A. V., & Claria, J. J. 2007, in Stellar 
Populations as Building Blocks of Galaxies, ed. A. Vazdekis, & R. F. 
Peletier (Cambridge: Cambridge University Press), 256 

Reddy, B. E., Tomkin, J., Lambert, D. L., & Allende Prieto, C. 2003, MNRAS. 
340, 304 

Rieke, G. H. & Lebofsky, M. J. 1985, ApJ, 228, 618 
Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 
Skrutskie, M.F., et al. 1997, in: F. Garzon et al., (eds.), The Impact of Large Scale 
Near-IR Skv Surveys, Kluwer Academic Publishing Co., Dordrecht, pp. 25. 
Sneden, C. 1973, ApJ, 184, 839 

Sulentic. J. W. & Tifft, W. G. 1973, The Revised New General Catalogue 
of Nonstellar Astronomical Objects (Tucson, Arizona: The University of 
Arizona Press) 

Tonry, J. & Davis, M. 1979, AJ, 84, 1511 

Villanova, S., Carraro, G, de la Fuente Marcos, R., & Stagni, R. 2004a, A&A, 
428, 67 



Villanova, S., Carraro, G, de la Fuente Marcos, R., 2004b, Proceedings of ASP 

Conference #317, p.196 
Unsold, A. 1955, Physik der Sternatmospharen, 2nd ed., Springer, Berlin, 

Heibelberg, New York 
Wielen, R. 1975, in Dynamics of Stellar Systems, ed. A. Hayli (Dordrecht: 

Reidel), 97 

Yi, S. K., Kim, Y.-C, & Demarque, P. 2003, ApJS, 144, 259 

Zacharias, N., Urban, S. E., Zacharias, M. I., Wycoff, G. L., Hall, D. M., Monet, 

D. G, & Raffeity, T. J. 2004, AJ, 127, 3043 
Zacharias, N., Finch, C, Girard, T., et al. 2009, AJ, in press 



